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Using the explicit form of the known single distribution functions (the Green's functions) in the 
double logarithm approximation of perturbative QCD, we analyze the structure of splitting diagrams 
as a source of double parton perturbative correlations in the proton. The related phenomenological 
effects are discussed for the conditions of the LHC experiments. 
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I. INTRODUCTION 

Strong interest has arisen [l| in the investigations and 
measurements of the multiparton interactions in high en- 
ergy hadron-hadron collisions. The analysis of final states 
with four jets, 7 + 3 jets, and W + 2 jets, performed by 
the AFS 0, UA2 Q, CDF 0,11, DO 0, and ATLAS 
Collaborations, provides convincing evidence for the sig- 
nificance of hard multiple parton interactions in these 
collisions and thereby supplements our understanding of 
the proton structure with new information. Studies of 
hard double parton scattering (DPS) have a long his- 
tory theoretically, with many references to prior work 
listed, for instance, in the recent review [l|. A greater 
rate of events containing multiple hard interactions is 
anticipated at the LHC, with respect to the experiments 
mentioned above, due to the much higher luminosity and 
greater energy of the LHC. Moreover, the DPS processes 
can constitute an important background [8|, |9| to signals 
from the Higgs and other interesting processes. Besides, 
certain types of multiple interactions will have distinctive 
signatures pi]|4l4l | , facilitating a detailed investigation of 
these processes experimentally and revealing information 
about parton pair correlations in the proton. 

The inclusive cross section of a DPS process in a 
hadron collision with the two hard parton subprocesses 
A and B may be written in the factorized form [see, for 
instance, Ref. [l3, where Eq. ([T]) is derived in detail] as 
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Here rij(a;i, X2;q; Qi, Q2) ^i'^ the generalized double 
parton distribution functions, depending on the longitu- 
dinal momentum fractions xi and X2 of the two partons 
undergoing the hard processes A and B at the scales 
Qi and Q2- ct,^. and o-^ are the parton- level subpro- 
cess cross sections. The factor m/2 is a consequence of 
the symmetry of the expression for interchanging parton 



species i and j. m = 1 li A = B and m = 2 otherwise. 
Note that these distribution functions also depend on the 
transverse vector q which is equal to the difference of the 
momenta of partons from the wave function of the collid- 
ing hadrons in the amplitude and the amplitude conju- 
gated. Such dependence arises because the difference of 
parton transverse momenta within the parton pair is not 
conserved. This transverse momentum q is the Fourier 
conjugated variable of the parton pair transverse separa- 
tion. The starting cross section formula ([T]) is somewhat 
similar to that usually used for single parton scattering. 
It was found (derived) in many works using the light-cone 
variables and the same approximations as those applied 
to the processes with a single hard scattering. 

The main problem is to make the correct calculation 
of the two-parton functions F^ (xi,a;2; q; Q^, (32) with- 
out p^l - l20| additional simplifying assumptions (the fac- 
torization of the impact-parameter dependence and x de- 
pendence that, by no means, should be treated as invi- 
olate). These functions and the corresponding evolution 
equations were considered in the current literature |2l| - 
[24| only for q = in the collinear approximation. In this 
approximation the two-parton distribution functions. 
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with the two hard scales set equal, satisfy the generalized 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
evolution equations, derived initially in Refs. [2l|, [23 |. 
Likewise, the single distributions satisfy more widely 
known and often cited DGLAP equations [25l - [28| . The 
functions in question have a specific interpretation in the 
leading logarithm approximation of perturbative QCD: 
they are the inclusive probabilities that in a hadron h 
one finds two bare partons of types i and j with the 
given longitudinal momentum fractions xi and X2- 

Based on these well-known collinear distributions, we 
have recently suggested [1^ a practical method which 
makes it possible to estimate the inclusive cross section 
for a DPS process without the oversimplified additional 
factorization assumption for rij{xi,X2; q = 0; Qi, Q2) = 
Dl^{xi',Qi)Dj^{x2',Q2) (which, in general, is inconsistent 
with the QCD evolution) but taking into account the 



QCD evolution explicitly [43|. Afterwards, similar re- 
sults were obtained also in Ref. |29j . with an emphasis 
on the differential cross sections, and were partly sup- 
ported in Ref. [30| , albeit with some diversity of opinion 
regarding a terminology mainly. We found that single 
and double perturbative splitting graphs can meaning- 
fully contribute to the inclusive cross section for a DPS 
process, in comparison with a "traditional" factorization 
component. 

The main purpose of the present paper is to analyti- 
cally study the structure of these single and double per- 
turbative splitting diagrams as a source of parton pair 
perturbative correlations in the proton. The paper is or- 
ganized as follows. In order to be clear and to introduce 
the denotations, we briefly recall some basic formulas 
from our previous work [l6| in Sec. II. The double par- 
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ton correlations in the double logarithm approximation 
of perturbative QCD arc estimated in Sec. III. The pos- 
sible phenomenological issues at the LHC are discussed 
in Sec. IV, together with conclusions. 



II. INCLUSIVE CROSS SECTION IN TERMS 
OF COLLINEAR DISTRIBUTIONS 



The inclusive cross section for DPS can be presented 
in the following form [16|: 
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or in substantially shorter yet less transparent form. 
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and for the combined ("interference") contribution. 
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Here asik'^) is the QCD coupling, and D^\{z;k'^,Q'^) 
are the known single distribution functions (the Green's 
functions generated by the usual DGLAP kernels) at the 
parton level with the specific 5-like initial conditions at 
Q^ = fc^. The one-parton distribution (before split- 
ting into the two branches at some scale fc^) is given by 
Dj^ {zi + Z2; /i^, fc^). The splitting functions 
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arc the nonrcgularized one-loop well-known DGLAP ker- 
nels without the "-I-" prescription. The single parton 
distribution functions D\^{xi] iJ.^,Q^) are the solutions of 
the DGLAP equations with the given initial conditions 
D\ {xi ; /x^) at the reference scale /x^ and may be expressed 
via the Green's functions Dl,{z;k'^,Q^) in the following 
way: 

1 
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Let us first consider the 1x1 component which de- 
scribes the two hard subprocesses A and B caused by the 
interactions of two pairs of partons in two independent 
branches of parton cascades. The probability of this dou- 
ble parton interaction depends on the spatial distribution 
(in the impact-parameter transverse plane) of these two 
branches of the parton cascade. In the momentum rep- 
resentation the spatial distribution is regulated by the 
two-parton (in a low x region this is mainly the two- 
gluon) form factor F2g{q). After the integration over q it 
gives the factor 



p4 f N_^ 



1 

feff 



which characterizes the transverse area occupied by the 
partons participating in hard collision and is often de- 
noted as an effective cross section, (Toff- 

Thus the value of the DPS cross section depends on 
the spatial correlations between the two partons in the 
incoming proton (hadron) wave function. Because of a 
strong k± ordering during the DGLAP evolution, the po- 
sition of the parton in the impact-parameter, bt, plane is 
frozen, and the form factor F2g describes the initial bt 
distribution formed in the nonperturbative region some- 
where at a low scale, less than jj?, where the DGLAP 
evolution starts. 

However, there is another sort of correlation caused 
by the splitting of one branch of the parton cascade into 
two branches. A splitting at the scale k^ produces two 
branches placed rather close to each other, with the spa- 
tial separation Sb^ ^ l/k'^. This effect is of perturbative 
origin and it may noticeably enlarge the DPS cross sec- 
tion, leading to a lower mean value of a^s- 



Depending on the kinematics of the DPS process and 
experimental cuts, we may concentrate on the two-parton 
correlations coming from the nonperturbative region, 
that is, on the 1x1 contribution or on the correlations 
of perturbative origin, like that in the 2x2 term. 

The contribution of the combined component (1 x 2) is 
regulated by the form factor i^2g (<?) [i3, ^M from the side 
of one incoming proton and by the perturbative splitting 
on another side. Since the form factor F2g{q) already 
provides the convergence of the q integral in the low q^ < 
/i^ region, here we mainly deal with the "long distance" 
correlations from the nonperturbative region. 

The nonperturbative correlations corresponding to a 
low q^ were discussed in detail in Refs. [13, 123, where 
it was proposed to consider the DPS events with small 
transverse momenta of the systems A and B produced by 
the hard subprocesses. This would be a very interesting 
study, albeit such cuts select a very small part of the 
total DPS cross section. An alternative possibility to 
study the correlations at a low scale q^ is to measure the 
asymmetric DPS processes where one "hard" scale Q\ is 
relatively low, say, the Xc-nieson and hi^-Ex dijet DPS 
processes. Since q^ < Q\ (Qi is a relatively low scale 
corresponding to Xc production) we have practically no 
space for pQCD splitting, and the 1x1 configuration will 
dominate. 

There is some discussion in the current literature [l9J, 
[2^, [so, [331 concerning the 2x2 component. This contri- 
bution incorporates the two splitting functions and the 
integration over q without the strong suppression factor 
F2g{q). Formally, in the region of not too small x, within 
the coUinear approach this contribution should be con- 
sidered as a result of the interaction of one pair of partons 
with the 2 ^ 4 hard subprocess [H, [H, [sO, [111 , since the 
dominant contribution to the phase space integral comes 
from a large q^ ^ min(Qf , Q|). However, as we argue in 
Ref. [i^ there may be configurations with a rather large 
interval between the splitting point (/c^, zi + Z2) and the 
(momentum) coordinates (Qf, xi) of the hard subpro- 
cess. In such a case the interval will be filled by the evolu- 
tion [either DGLAP when fc^ « Ql or Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) ^^^ when Zi + Z2 » Xi] 
which will produce additional secondaries. The corre- 
sponding process can not be described by the 2^4 
hard matrix element. Here we have to use our fornui- 
las @ and Q, especially in the case of a configuration 
with two quite different scales (for instance, Qf << Q\)- 

In order to understand the structure of these additional 
contributions (1x2 and 2x2) better, we consider them 
in the double logarithm approximation, which allows us 
to obtain some analytical estimations. 



III. DOUBLE LOGARITHM APPROXIMATION 

Let us write down all the integrations with splitting 
functions separately to make the analysis more transpar- 



ent, 
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Because of the strong suppression faetor F2„{q) in a 
single splitting diagram (1x2 contribution), this inte- 
gral Z?^^2(^i'2;2;9^,Qi, Q2) can be estimated at the ref- 
erence scale q^ = fj,^ and can be considered as the QCD 
evolution correction to the factorized double parton dis- 
tribution functions. For double splitting diagrams (2x2 
contribution) we should keep the further nonlogarithmic 
integration over q in mind. 

In the double logarithm approximation we can re- 
strict ourselves to the main gluon contribution only and 
rewrite the integral under consideration in the following 
form Q: 

D^^2{x,,X2-,q\QlQl) 
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where u = zi + Z2 and z = zi/u. 

The limits in the u and z integrations are Xi < uz, 
X2 < u(\ — z), u < 1, and z < 1. The Green's functions 
(gluon distributions at the parton level) in the double 
logarithm approximation (see, for instance, J27ll33l|) read 



(10) 



xD'^g{x, t) ~ ANctv'^^"^ exp [v - at]/\/27r, 
where v = yJ%Nct\n{l/x), a = ^-N^ + ^rif/N^, 
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and where (3 ~ (lliVc— 2n^)/3, with the number of active 
flavors Uf, A is the dimensional QCD parameter, and 
iVc = 3 is the color number. Recall that in Eq. (|TT|) the 
one-loop running QCD coupling 
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/31n(g2/A2) 
was used. After that, the integral ^ may be rewritten 
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where 



d{t,y,z)= Vty+ViTi-t){Yi- 
+ ^{T2-t){Y2-v) 
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with t = t{P), Ti = t{Ql), T2 = t{Ql), T = f(g2); and 
y = ln(l/u), Yi = ln(l/a;i) - ln(l/z), Y-, = ln(l/a;,) - 
ln(l/(l-z)). 

In Eq. (J13p we keep the leading exponential terms only, 
which have the same structure both at the parton level 
and at the hadron level under smooth enough initial con- 
ditions at the reference scale. Indeed, in the double log- 
arithm approximation Eq. ([7]) reads 
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with T = ^(Q^) and Y = ln(l/x). The y' integration is 
not a saddle-point type, and therefore, one of the edges, 
namely y' — )■ {z' — >■ 1), dominates, provided that the 
initial gluon distribution does not grow too much as z' de- 
creases. In fact, one needs z'Df^{z', 0) ^ {l/z')°- at z' ^- 
with a < A, where A = yj2NcT/Y > 0. Note that the 
parametrization of the initial gluon distributions usually 
used satisfies this condition (e.g., the CTEQ parametriza- 
tion from Ref. j37|). 

We are interested in the domain with large enough Ti , 
T2, ln(l/xi), and ln(l/x2), when the exponential factors 
are large in comparison with 1 and where the approxima- 
tions above are justified. In this case the integration over 
the rapidity y has a saddle-point structure in the wide in- 
terval of z integration not near the kinematic boundaries. 
The saddle-point equation reads 



Vt v/(Ti - t) v/(T2 - t) 
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It may be solved explicitly in the simplest case of the 
two hard scales set equal, Ti = T2 — T, and at Yi ~ 
Y2 ^Y ^ ln(l/.T), i.e., in the z region where ln(l/z) << 
ln(l/a;) and ln(l/(l - z)) << ln(l/x) [H. Then the 
saddle point is equal to 



yo = Yt/{4T ~ 3t) 



(17) 



and Eq. ([T3|) reduces to 
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The t integration is not a saddle-point type, and there- 
fore, one of the edges, namely t ^ t, dominates. That 

is. 



x^Dll{x,x;T,T,T)r^cxp[^/mc^/Y(4T^^3?j]. (19) 



IV. DISCUSSION AND CONCLUSIONS 

Now let us discuss in more detail what follows from our 
estimation of splitting integrals in the double logarithm 
approximation by the saddle-point method [46[. 

For single splitting diagrams (1x2 contribution) the 
lower limit for the t integration in the estimation (|18p 
may be taken at the reference scale, i.e., r = t((7^)|g=^ = 
0, due to the strong suppression factor F2g{q). The char- 
acteristic value of q is of the order of the "effective gluon 
mass" nig ^ /i in the further q integration. Thus, one 
obtains for this contribution the following estimation: 



'Dllix,x;0,T,T)^exp[ 



'Tt)]. (20) 



This means that the splitting takes place at the "charac- 
teristic point" with the scale fc^ close to fi^ and with the 
longitudinal momentum fraction u ^ 1 (the saddle point 
ya ^ t ^ T ^ in this case). After splitting, one has two 
independent ladders with the well-developed BFKL [33- 
i^ and DGLAP [H-il evolutions. Every ladder con- 
tributes to the cross section with the large exponential 
factor cxp [-i/SN^VYT] , which is just the same as for 
single distributions [compare Eq. (pil)) with Eq. pH]) ]. 
Therefore, in the double logarithm approximation sin- 
gle splitting diagrams © have, in fact, the factorization 
property if one takes only the leading exponential fac- 
tors into consideration. However, the contributions to 
the cross section from single splitting diagrams and from 
the factorization component differ in nonexponential fac- 
tors omitted here, especially, in the different "normal- 
ization" at the reference scale: one "nonperturbative" 
parton for single splitting diagrams and two initial inde- 
pendent "nonperturbative" partons for the factorization 
component. 

The factorization property for the integral 
D]^2{xi,X2]fJ-^,Qi,Q2) (0) was found also in Ref. [Sg in 
the double logarithm limit based on other techniques. 

For double splitting (2 x 2) diagrams the leading expo- 
nential contribution arises from the lower limits of t inte- 
gration and either lower or upper limits of q integrations, 
depending on the available rapidity interval Y. There is 
competition between the exponential factor caused by the 
evolution, which prefers a small r, and the phase space 
factor in the q^ integral. Because of the nonlogarithmic 



character of the integration over (Pq for a not sufRciently 
large Y, the contribution from the upper limit of q may 
dominate. Indeed, let us consider the production of two 
bb pairs in a central rapidity (rj ^ 0) region. That is, we 
take Ti = T2 = T,Yi ~ Y2 = Y and keep just the leading 
exponential factors in the double parton distributions, 

x^Dh2{x,x,q'^,Q'^,Q'^) 

~ exp(V8iVc>"(4T - 3t) - 2aT + ar) , (21) 

where for better accuracy we keep the term —at in the 
exponent of Eq. (fTU]) [recall that here t{q^) = r]. 
Thus the logarithmic dq^ /q'^ integral takes the form 

/ ^ exp(2 v'8A^c^(4r - Sr) - AaT + 2aT)q^ , (22) 

with ln(gVA2) = L = \n{n^ /A^)e'^''/^. The L behavior 
of the integrand of Eq. ([22]) . 
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is shown in Fig. 1 in the case of F = 5 and Y — 2, 
corresponding to the LHC energy -y/s — 14 TcV and the 
RHIC energy y/s = 500 GeV (to be more or less realistic, 
here Y is calculated a.sY = ln(xo/a;) with xq ^ 0.2). For 
this numerical estimation we take Nc = 3, 71/ = 4 and 
A = 150 MeV [13; Q2 = 50 GeV^. In the LHC case 
we consider also the DPS ly-boson production, taking 
g2 = 10* GeV^ and Y = 3. 

As it is seen in Fig. 1, where the relevant quantity 
a1{q^) f {L) is plotted, for the DPS production of two hh 
pairs the major contribution comes from a low q^ . That 
is, the reaction may be effectively described by the 1x1 
term; the formation of two parton branches (one to two 
splitting) takes place mainly at low scales. However, at 
the RHIC energy, when the available rapidity interval is 
not large, the q^ dependence is not steep and the contri- 
bution caused by the splitting somewhere in the middle of 
the evolution is still not negligible. The same can be said 
about the DPS M^-boson production at the LHC. Here 
the upper edge of the q^ integral dominates. This part 
may be described as the collision of one pair of partons 
supplemented by a more complicated, 2 — >• 4 or 2 — >• 2W , 
hard matrix element. However, clearly we need to also 
account for contributions from the whole q^ interval. 

In other words, depending on the precise kinematics, 
we may deal either with a single parton pair collision 
(times the 2 -> 4 hard subprocess), with the contribu- 
tion of the 1x1 type, where the formation of two par- 
ton branches (one to two splitting) takes place at low 
scales, or with the 2x2 configuration where the splitting 
may happen everywhere (with more or less equal prob- 
abilities) during the evolution. Note that just this last 
possibility may be relevant for the LHC experiments. 
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"W LHC 



bB RHIC 



q^ (GeV^) 

FIG. 1: The q dependence of the integrand f{L) in the 
logarithmic scale. 



with jets, four-jet production, production of a hh pair 
with two jets, production of a 65_pair with W , have been 
discussed in detail in a review [ij, with many references 
to prior works listed therein. Quite recently the processes 
with pairs of heavy quarkonia in the final state were con- 
sidered [39l - l4]| as precise probes of the DPS at the LHC. 
We believe that production of a hh pair (or J/V') with 
W may be a good candidate process to probe the QCD 
evolution of the double distribution functions due to a 
configuration with two quite different scales. 

For completeness, it is also interesting to estimate the 
value of the exponential factors available at the LHC 
kinematics in single ladder diagrams. The asymptotic 
behavior of the distribution functions is determined by 
the factor 



exp [^WcWt] = exp [2.4v/fcBFKLv/fcDGLAp] , (24) 
where 



fcBFKL = l^ = ln(l/x)<14 



cLt Xti 
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When both splittings take place at relatively small 
scales ^ fi^ << Q^, the 2x2 contribution cannot be 
considered as the result of the interaction of one pair 
of partons with the 2 — )■ 4 hard subprocess, unlike the 
statement in Rcfs. 0, [H, HO, [ll- 

Recall also that there may be a configuration with two 
quite different scales (say, Qf << Q2), in which the up- 
per limit of the q^ integral is given by a smaller scale 
{a.t q > Qi the hard matrix element corresponding to a^ 
begins to diminish with qt). In this case the collinear 
evolution from the scale q = Qi up to the scale Q2 
in the ladders is sufficiently justified. A configuration 
with two quite different "final" rapidities [for instance, 
1 << ln(l/a;i) << ln(l/x2)] is also interesting to probe 
single and double splitting diagrams, since the BFKL 
evolution takes place in the ladders before and after (in 
one of ladders only) splitting. 

Here it is worth noticing that the asymptotic predic- 
tion mainly "teaches" us a tendency and tells us nothing 
practical about values of xi, X2, Qi, and Q2, beginning 
from which the asymptotic behavior is a good approx- 
imation to the real one. Therefore, it makes sense to 
consider the quantitative contribution of the 2x2 term 
even within the collinear approach as applied to the LHC 
kinematics, where the large available values of Qi and Q2 
(in comparison with rUg and /z), ln(l/xi) and ln(l/x2) 
(in comparison with 1) can provide configurations with 
the BFKL or DGLAP evolution in ladders before and 
after splitting, depending on the processes under consid- 
eration. 

A number of processes were suggested in order to probe 
DPS at the LHC. Promising candidate processes, such 
as same-sign W production, Z production in association 



fcoGLAP = ln[(ln(Q/A)/(ln(M/A)] < 1.4 

at g,„ax = 250 GeV, ^i = 1 GcV, A = 0.15 GeV. In 
Eq. ^ we put uj = 4, A^c = 3. 

However, it is better to compare the pure BFKL factor 
LUo ln(l/a;) ~ {asNc/ir) ln(l/x) with the analogous 
DGLAP factor 2Nct. In the BFKL case the resum- 
mation of the next-to-leading logarithmic corrections 
leads to (jJq ^ 0.3 [42|, that is, the BFKL power 
luq * ln(l/a:;) < 4, while for the DGLAP evolution we 
have (4A''c//3)A:dglap < 2. Thus the LHC kinematics 
admits a wider interval for the BFKL evolution than for 
the DGLAP one in compatible dimensionless variables. 
Therefore, it may be interesting (and even more justified 
theoretically) to consider the multiple parton interac- 
tions in the framework of the BFKL approach (sec. for 
instance, Ref. J2G|). 

In summary, we have demonstrated that all compo- 
nents of the generalized double distribution functions 
have the factorization structure in the double logarithm 
approximation and contribute to the cross section with 
the same leading exponential terms in Y and T, but with 
different weights (nonexponential factors). For the de- 
batable double splitting diagrams, depending on the pre- 
cise kinematics, we may deal either with a single parton 
pair collision (times the 2-^4 hard subprocess) , with the 
contribution of the 1x1 type where the formation of two 
parton branches (one to two splitting) takes place at low 
scales, or with the 2x2 configuration where the splitting 
may happen everywhere (with more or less equal proba- 
bilities) during the evolution. In order to probe the QCD 
evolution of the double distribution functions better, we 
suggest to also investigate the processes with two quite 



different scales, in particular, production of a bb pair (or 
J/tp) with W, which was estimated at the LHC kinemati- 
cal conditions in Ref. [IJI using the factorized component 
only. 
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